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Abstract 

The cross section for exclusive heavy quark and heavy antiquark pair {QQ) production in periph- 
eral ultrarelativistic heavy ion collisions is calculated for the LHC energy ^/snn = 5.5 TeV. Here 
we consider only processes with photon-photon interactions and omit diffractive contributions. We 
present results in the impact parameter equivalent photon approximation (EPA) and compare some 
of them with results obtained by exact calculations of the Feynman diagrams in the momentum 
space. We include both QQ, QQg and QQqq final states as well as photon single-resolved compo- 
nents. Realistic charge densities in nuclei were taken in the calculation. The different components 
give contributions of the same order of magnitude to the nuclear cross section. The cross sections 
found here are smaller than those for the diffractive photon-pomeron mechanism and larger than 
diffractive pomeron-pomeron discussed in the literature. 
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I. INTRODUCTION 



Heavy quark - heavy antiquark production was studied in the past in photon-photon, 
photon-proton, inclusive and exclusive proton-proton and heavy ion collisions when the ions 
break apart. In principle, the heavy quark - heavy antiquark pairs can be produced also 
in exclusive coherent 77 processes when the nuclei stay intact. These processes, in analogy 
to lepton pair production |2[, should be "increased" by the large charges of the colliding 
nuclei. On the other hand the effects of the nuclear form factors which diminish contribu- 
tions of large energy virtual photons should be included. Recently two of us (M.K.+A.S.) 
have shown that the inclusion of realistic form factors (corresponding to realistic charge 
densities) is essential for reliable estimation of the cross sections for exclusive lepton-pair 
production especially for large rapidities and large transverse momenta of the produced par- 
ticles. The inclusion of the realistic form factor should be even more important for heavy 
quark - heavy antiquark pair production. 

In contrast to dileptons the QQ state cannot be directly observed. In practice one mea- 
sures rather heavy mesons or electrons from their semileptonic decays. Then the final states 
are already complicated due to hadronization process. Therefore one has to include also 
different, yet simple, partonic states as QQg and QQqq. 

Recently one of us (M.M.) with coworkers made an estimation of the cross sections for 
diffractive mechanisms [J]. It is of interest to make a realistic estimate of the cross section 
for photon-photon mechanism and compare it with the diffractive contribution. 

II. HEAVY QUARK HEAVY ANTIQUARK PAIR PRODUCTION 
A. Photon photon subprocesses 
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FIG. 1: Representative diagram for the Born amplitude. 

In the present analysis we include the mechanisms shown in Fig. [H [2l |3l HI The first 
one (called sometimes direct) is identical as the one for the production of charged lepton 
pairs. In contrast to the dilepton production in the case of the quark production one has to 
include also QCD corrections. Corresponding diagrams for NLO approximation are shown 
in Fig. |2j Heavy quarks can be also produced in association with light quark- antiquark pairs 
as shown in Fig. |3l The last two diagrams correspond to single-resolved components when 
only small part of one photon interacts with the other photon. All these processes were 
studied in detail in Refs.j5|-0]. 

Let us start with the Born direct contribution. The leading-order elementary cross section 
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FIG. 2: Representative diagrams for the leading-order QCD corrections. 
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FIG. 3: Representative diagrams for QQqq production. 



for 77 — i- at 2-photon energy W^^ takes the simple form 
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where QQ = cc, bb, A^c = 3 is the number of quark colors, w = y 1 — -p^T^ and eg is the 

fractional charge of the heavy quark. The formula has been derived for the first time in Ref. 
[sl. In the current calculation we take the following heavy quark masses: rric = 1.5 GeV, 
rub = 4.75 GeV. This formula can be directly used in the impact-parameter-space (called 
here fe-space for brevity) equivalent photon approximation (EPA), as we will see below. It 
is obvious that the final QQ state cannot be observed experimentally due to the quark 
confinement and rather heavy mesons have to be observed instead. It is noticed that the 
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FIG. 4: Representative diagrams for the single-resolved mechanism. 



presence of additional few light mesons is rather natural. ^ This forces one to include more 
complicated final states. 

In contrast to QED production of lepton pairs in photon-photon collisions, in the case of 
QQ production one needs to include also higher-order QCD processes which are known to 
be rather significant. Here we include leading-order corrections only for the dominant, in 
heavy-ion collisions, direct contribution. The details concerning the higher-order corrections 
to heavy quark and heavy antiquark production in photon-photon collisions can be found 
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In a^-order there are one-gluon bremsstrahlung diagrams (77 — t- QQg) and 
interferences of the Born diagram with self-energy diagrams (in 77 — )■ QQ) and vertex- 
correction diagrams (in 77 — )■ QQ). The relevant diagrams are shown in Figj2l In the 
present analysis we follow the approach presented in Ref. 
be written as ^ 
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(2.2) 



The function f^^^ is calculated using a code provided by the authors of Ref. [IT!] which uses 
program package HPL In the present analysis the scale of a.^ is fixed at = 4mQ. 

We include also the subprocess 77 — )■ QQqq, where q (g) are light, u, d, s, quarks (anti- 
quarks). The cross section for this mechanism can be easily calculated in the color dipole 
framework [SirH]- In the dipole-dipole approach j6| the total cross section for the 77 — )■ QQ 
production can be expressed as 
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where (p, z) are the quark - antiquark wave functions of the photon in the mixed 
representation and add is the dipole-dipole cross section. Eq. (12.31) is correct at sufficiently 
high energy W^y^ ^ 2mQ. At lower energies, the proximity of the kinematical threshold is a 
concern. In Ref. ^ a phenomenological saturation-model inspired parametrization for the 
azimuthal angle averaged dipole-dipole cross section has been proposed: 
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Here, the saturation radius is defined as 
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and the parameter Xab which controls the energy dependence was given by 
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^ There are also final states with exclusively two Qq and Qq mesons 
are however much smaller. 



The corresponding cross sections 
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In the numerical calculations we are using the model parameters (xq, A, cxo 
effective radius r^g = {P1P2Y / {pi + P2) • Some other parametrizations of the dipole- 
dipole cross section were discussed in the literature (see e.g. [loj). The cross section for 
the 77 — )■ QQqq process here is much bigger than the one corresponding to the tree-level 
Feynman diagram 15| as it effectively resums higher-order QCD contributions. 

As discussed in Ref. ^ the QQqq component have very small overlap with the single- 
resolved component because of quite different final state, so adding them together does not 
lead to double counting. The cross section for the single-resolved contribution can be written 
as: 
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and a^g are ele- 
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where gi and g2 are gluon distributions in photon 1 or photon 2 and 
mentary cross sections. In our evaluation we take the gluon distributions from Ref. 

In Fig. [5] we show the elementary cross sections for all processes as a function of the 
photon-photon center-of-mass energy. For charmed quark the direct term dominates at 
low energies near the threshold while the four-quark component at slightly larger energies 
and the resolved components at even larger energies. For bottom quarks the four-quark 
component is always larger than the direct term. 
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FIG. 5: The elementary cross section for the different processes as a function of the photon- 
photon center-of-mass energy. In the left panel we show the Born cross section (solid line) and 
leading-order QCD corrections (dash-dotted line). In the middle panel we show contribution of 
four-quark final states as calculated in the saturation model. In the right panel we show in addition 
contribution of single-resolved processes. 



It is not clear a priori how this will change in the nucleus-nucleus collisions where one 
should take into account photon-photon luminosities. This will be discussed in the result 
section. 



B. 6-space EPA 



Here we wish to sketch the 6-space EPA used in the present analysis. The details on its 
development can be found in Ref. |3|]. 
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The total cross section for the PbPb — )■ PbPbQQ process can be factorized into the 
equivalent photon spectra, {u, b), and the 77 — )■ QQ subprocess cross section as: 

a (PbPb ^ PbPbQQ- spbPfe) = la{-f-f^ QQ; W^^) 9 (|bi - bal - 2Ra) 

X (wi, bi) N (^2, bs) d^bid^badwidwa • (2.8) 

After performing a change of integration variables, the cross section can be expressed as the 
five-fold integral: 

a (PbPb ^ PbPbQQ; spbPfc) = 1 a (77 ^ QQ\ W^^) 9 (|bi - bsl - 2Ra) 

xN {uji,hi)N {u2,h2)27ibdbdb^dby^^dW^^dY , (2.9) 

where the quantities b^ = (bi^ + b2x)/2, by = {biy + &2s/)/2 and b = bi — b2 have been 
introduced. Eq. (12. 9 p is used to calculate the total cross section for the PbPb — )■ PbPbQQ 
reaction as well as the distributions in impact parameter b = |b|, W^^ = Mqq and quark 

pair rapidity Y{QQ) = | (i/q + uq^. A detailed derivation of formula f l2.9p can be found in 

[3|. The photon flux can be expressed in terms of the charge form factors F{Q'^) as: 



where Ji is the Bessel function of the first kind, Q^ = —q^ > and q is the four-momentum 
of the quasireal photon. The form factor is the Fourier transform of the nucleus charge 
distribution, p(r): 

Fm= /^p(r)sin(Qr)rrfr = l-^^ + ^^ + ... . (2.11) 



III. RESULTS 



Now we will discuss the nuclear cross sections obtained within 6-space EPA described 
shortly above. 

Let us start with the direct contribution. In Fig. Elwe show the differential distribution in 
photon-photon energy which for the direct component is also the distribution in the quark 
- antiquark invariant mass. They are presented without (dashed line) and with (solid line) 
absorption effects. This is done in the 6-space by integrating over impact parameter either 
in the full range of 6 > (without absorption) or for b > Ri + R2 (with absorption), where 
Ri^2 are the nuclei radii. For lead nucleus Ra = 1.2 A^^^ ~ 7 fm. The results for the cc 
are shown on the left hand side whereas the results for the bb on the right hand side of the 
figure. One can clearly see that the cross section for the cc pair is considerably larger than 
that for the bb pair (this is because of the charge and mass differences). We will return to 
it when discussing total cross sections at the end of this section. 

Above we have shown results obtained in the equivalent photon approximation. The same 
observable can be also obtained calculating Feynman diagrams in the momentum space. The 
details of this calculation have been carefully presented in paper jsf. In Fig. [7] we compare 
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FIG. 6: The 77 subsystem energy distribution, , for PbPb — )• PbPbcc (left panel) and PbPb — >■ 
PbPbbb (right panel). The solid line denotes the cross section calculated within EPA approach for 
peripheral collisions (6 > 14 fm) while the dashed line includes also central collisions. 
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FIG. 7: The nuclear cross section as a function of 77 subsystem energy for the PbPb — ?• PbPbcc (left 
panel) and for the PbPb — )• PbPbbb (right panel) reactions calculated in the EPA approximation 
(dashed lines) and in the momentum space (solid line). 



the EPA results and those for momentum space calculation. The numerical results for the 
two methods are quite similar. 

In Fig. [H] we compare the contributions of different mechanisms discussed in the present 
paper as a function of the photon-photon subsystem energy. For the Born case it is identical 
as a distribution in quark-antiquark invariant mass. In the other cases the photon-photon 
subsystem energy is clearly different than the QQ invariant mass. Therefore, this distribution 
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is rather theoretical and would be difficult to measure experimentally. These distributions 
reflect the energy dependence of the elementary cross sections (see Figj5]). Please note a 
sizable contribution of the leading-order corrections close to the threshold and at large 
energies for the cc case. Since in this case W^^ > Mqq, it becomes clear that the QQqq 
contributions must have much steeper dependence on the QQ invariant mass than the direct 
one which means that large QQ invariant masses are produced mostly in the direct process. 
In contrast, small invariant masses (close to the threshold) are populated dominantly by 
the four-quark contribution. Therefore, measuring the invariant mass distribution one can 
disentangle the different mechanisms. As far as this is clear for the cc it is less transparent 
and more complicated for the bb production. In the last case, the experimental decomposition 
may be in practice not possible. 




FIG. 8: The nuclear cross section as a function of photon-photon subsystem energy W^^ in EPA. 
The red solid line denotes the results corresponding to the Born amplitude (cc -left panel and bb 
-right panel). The leading-order QCD corrections are shown by the red (on line) dash-dotted line. 
For comparison we show the differential distributions in the case when an additional pair of light 
quarks is produced in the final state (dashed lines) and for the single-resolved components (dotted 
line). 

Another distribution which can be calculated in the 6-space EPA is the distribution in 
rapidity of the particles produced in the flnal state: i.e. rapidity of the QQ pair for the 
direct component, the rapidity of the whole QQg or QQqq system for the QQg and four- 
quark component or the rapidity of the very complicated system for the single-resolved 
components. In Fig. [9] we present distributions in such a variable for all components, for 
cc (left panel) and bb (right panel) production. It may be quite difficult to reconstruct 
the QQg and four-quark rapidity distribution experimentally and even more difficult to 
reconstruct the rapidity of the complex flnal state for the single-resolved components. The 
distribution for the single-resolved components is narrower than for the other components 
which is related to the fact that large part of the energy is taken by the photon remnants. 

Finally, in Tables 1, 2, 3 and 4 we have summarized the results for the total cross sections 
for different components calculated within distinct methods: 6-space EPA, exact momentum 
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FIG. 9: The rapidity distribution, in the 6-space EPA. Here the impact parameter is in the 
whole range of impact parameter (0 < 6 < cxd). The red (on hne) sohd hne denotes the results 
corresponding to the Born amplitude (cc -left panel and bb -right panel) and the dash-dotted 
line corresponds to the leading-order QCD corrections. For comparison we show the differential 
distributions in the case when an additional pair of light quarks is produced in the final state 
(dashed lines) and for the single-resolved component (dotted line). 

space and momentum-space EPA described in detail in j^. 

In Table 1 we present the cross sections for the Born direct component only. In order 
to illustrate the absorption effect we show both the integral calculated from to "infinity" 
and the integral calculated from Ri + R2 to "infinity" . The integration in b is only slowly 
convergent, especially for the lighter cc pairs. Therefore we also show the practical upper 
limit dependence of the cross section. We are not sure that the upper limit for cc is sufficient. 
The results obtained within momentum-space EPA j2| are very similar to those obtained in 
the exact momentum space method {3| . For bb these results are also in good agreement with 
the 6-space EPA. 

In Table 2 we have ensambled leading-order QCD corrections. They constitute about 
one third of the Born contribution. 

In Table 3 we show results for the QQqq components obtained using 6-space EPA. The 
cross sections here are of the similar order of magnitude as those for the direct component. 
For charm, the cross section for QQqq is approximately equal to the direct contribution. On 
the other hand, for bottom the result is almost four times bigger than the direct component. 
This is due the dominance of four-quark component even near threshold as shown in Fig. 
5. That feature is already known from Ref. j6[, where the dipole-model was compared with 
LEP data for (j(e+e- ^ bbX). 

For completeness in Table 4 we present results for the single-resolved component. The 
cross sections for this component is comparable to the other components. 

Finally, in Table 5 we present cross sections which include all discussed mechanisms and 
their relative contributions. We see that the Born mechanism dominates for cc production 
but four-quark component for bb production. 



9 



TABLE I: Total cross section for the QQ component calculated within different methods. 



Process 


6-spac 

b>0 


e EPA 

b>14 fm 


Momentum- 
-space 


Momentum- 
-space EPA 


PbPb-^ PbPbcc 
PbPb^ PbPbbb 


1.18 mb (b<4000 fm) 
1.13 mb (b<1000 fm) 
2.53 fib 


1.05 mb (b<4000 fm) 
1.00 mb (b<1000 fm) 
2.05 fib (b<1000 fm) 


1.36 mb 
2.54 fib 


1.230 mb 
2.54 fib 



TABLE II: The leading-order QCD corrections to the total cross section within 6-space EPA. 



Process 


6-spa( 
b>0 


ze EPA 
b>14 fm 


PbPb PbPbcc 


0.41 mb 


0.36 mb 


PbPb PbPbbb 


1.00 fib 


0.83 fib 



The event rates should be large also after hadronization process. For example, the total 
cc and bb two-photon rates in peripheral PbPb collisions over a 10^ s run at LHC are 
N{cc) = 10.4 X 10^ and N{bb) = 4.6 x 10^ including Born, QCD-corrections, single-resolved 
and four-quark component. This is done using the 6-space EPA with absorption and taking 
a luminosity of Cp^pb = 4.2 x 10^^ cm~^s~^. The results in Tables 1 and 4 can be compared 
directly to previous studies which rely on QCD collinear factorization approach, for instance 
in Ref. [2l|. Notice that for charm the direct contribution produces similar cross sections 
compared to j2l| whereas for bottom it is a factor two bigger. The QQqq contribution was 
not included in |2l|] and in addition the single and double resolved processes were shown to 
be negligible for LHC energy compared to the direct one. As a short final comment, the 
cross sections presented in Tables 1 and 4 are larger than the estimations for the exclusive 
heavy-quark production in double-pomeron exchange (DPE) process (it was found in Ref. 
[i*] crp^p^{cc) ~ 4.2 fih and (7p™b(^^) ~ 0-2 A^b) whereas they are smaller than the quark-pair 
production in photon-pomeron processes. 

IV. CONCLUSIONS 

In the present paper we have concentrated on production of heavy quark - heavy an- 
tiquark pairs in coherent photon-photon subprocesses. A discussion on similar diffractive 
processes have been already presented in the literature and will be not repeated here. The 
photon-photon processes are dominant in the case of exclusive production of charged lepton 
pairs. In our calculations we have used realistic nuclear form factors calculated as Fourier 
transform of the realistic charge density of the nucleus known from the electron scattering 

TABLE III: Total cross section for the QQqq components within 6-space EPA. 



Process 


6-spa 
b>0 


ze EPA 
b>14 fm 


PbPb PbPbccqq 


0.82 mb 


0.67 mb 


PbPb PbPbbbqq 


9.40 fib 


6.98 fib 
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TABLE IV: Total cross section for the single-resolved components within 6-space EPA. 



Process 


6-spa( 

1)>() 


ze EPA 

l)>l-i f m 


PbPb PhPh c- 


0.52 nib 


0.39 mb 


PbPb^ PbPbbb 


1.51 nh 


0.97 nb 





CABLE V: Partia 


contributions of different mechanisms. 




0"tot 


Born 


QCD-corrections 


4-quark 


Single-resolved 


cc 


2.47 mb 


42.5 % 


14.6 % 


27.1 % 


15.8 % 


bb 


10.83 nb 


18.9 % 


7.7% 


64.5 % 


8.9 % 



off nuclei. Recently, this was shown to be crucial for reliable estimation of the exclusive 
lepton pair production. 

We have calculated cross sections for exclusive production of charm-anticharm and 
bottom-antibottom pairs, for the QQg and QQqq final states as well as for the single- 
resolved components in the high-energy peripheral lead-lead collisions for the LHC energy 
V^NN ~ ^-^ GeV. Large cross sections have been found in the case of charm quarks (an- 
tiquarks) production. In contrast to the exclusive dilepton production in the case of the 
heavy quark - heavy antiquark production large QCD corrections appear. Their fractional 
contribution strongly depends on the photon-photon subsystem energy. 

Two methods have been used to calculate the Born contribution: impact parameter 
equivalent photon approximation (6-space EPA) and the Feynman diagrammatic approach 
for the Born component called here momentum space approach. The 6-space EPA is an 
approximation but allows to include absorption effects in a simple way by limiting the range 
of integration over impact parameter. The direct contribution was calculated in both the 
6-space EPA and the Feynman graph approach while the leading-order QCD correction, 
four-quark component as well as single-resolved components only in the 6-space EPA. 

Wc have presented total (phase space integrated) cross sections as well as some selected 
differential distributions relatively easy to calculate in the 6-space EPA. The absorption 
effects turned out to be larger for bottom quarks (20 %) than for charm quarks (10 %). 
Since both methods lead to similar effects one can use the momentum space approach to 
calculate, or at least to estimate, different observables which are not straightforward in the 
6-space approach. 

We have found that the contributions of two- and four- quark and single-resolved final 
states are of similar size. We have found also that the large invariant masses of the QQ system 
are populated dominantly by the direct 77 — > QQ subprocesses while smaller invariant 
masses by the 77 — >■ QQg, 77 — )■ QQqq or single-resolved components. This could be 
potentially helpful in experimental identification of the all components. There are known 
experimental methods how to distinguish large transverse momentum b (6) jets, therefore 
exclusive measurement of such jets should be possible in the LHC experiments. 

The cross sections found for the QED processes discussed in the present paper seem 
smaller than those found in diffractive photon-pomeron processes but smaller than diffrac- 
tivc pomcron-pomcron process. The diffractive processes are more difficult to be reliably 
calculated. A combined simultaneous analysis of all processes including different differential 
distributions seems indispensable in the future. Since the QED processes, as demonstrated 



11 



here, can be reliably calculated it can be used as a background to the much more involved 
diffractive processes. 
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